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Abstract 

We study the robustness of the selftrapping phenomenon exhibited by the Discrete Non- 
linear Schrodinger (DNLS) equation against the effects of nonadiabaticity and quantum 
fluctuations in a two-site system (dimer). To test for nonadiabatic effects (in a semiclassi- 
cal context), we consider the dynamics of an electron (or excitation) in a dimer system and 
coupled to the vibrational degrees of freedom, modeled here as classical Einstein oscillators 
of mass M. For relaxed (coherent state) oscillators initial condition, the DNLS selftrap- 
ping transition persists for a wide range of M spanning 5 decades. When undisplaced 
initial conditions are used, the selftrapping transition is destroyed for masses greater than 
M ~ 0.02. To test for the effects of quantum fluctuations, we performed a first-principles 
numerical calculation for the fully quantum version of the above system: the two-site Hol- 
stein model. We computed the long-time averaged probability for finding the electron at the 
initial site as a function of the asymmetry and nonlinearity parameters, defined in terms of 
the electron-phonon coupling strength and the oscillator frequency. Substantial departures 
from the usual DNLS system are found: A complex landscape in asymmetry-nonlinearity 
phase space, which is crisscrossed with narrow "channels", where the average electronic 
probability on the initial site remains very close to 1/2, being substantially larger outside. 
In the adiabatic case, there are also low-probability "pockets" where the average electronic 
probability is substantially smaller than 1/2. 
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1 Introduction 



The Discrete Nonlinear Schrodinger (DNLS) equation is a paradigmatic equation that de- 
scribes properties of chemical, optical or condensed matter systems where selftrapping mech- 
anisms are present. It can describe the dynamics of a set of anharmonic oscillators Q, the 
power exchange between nonlinear optical couplers P] or the motion of a quantum mechan- 
ical particle propagating in a discrete medium while interacting strongly with vibrations ||. 
In this last case, for an electron (or excitation) propagating in a discrete one-dimensional 
medium, DNLS has the form 

i^ = V(C n+1 + C n „ 1 )- X \C n \ 2 C n (1) 

where C n (t) is the probability amplitude of finding the electron on site n at time t, V is 
the nearest-neighbors hopping term and x is the nonlinearity parameter, proportional to 
the square of the (strong) electron-phonon coupling on site n. The DNLS equation can be 
"derived" as the antiadiabatic limit of the problem of a quantum excitation propagating in 
a lattice and coupled to the vibrational degrees of freedom of each site, modelled as a set 
of classical Einstein oscillators 

dC 

i_ df = V ( Cn+1 + Cn ~^ + aU nCn (2) 

M u n + k u n = — a \C n \ 2 (3) 

After the transients, the oscillator amplitude u n (t) is given byu n = (l/y/Mk) J^dt' sin(a>(i- 
t') a\C n (t')\ 2 . If we assume the oscillator time scale much shorter than the electronic time 
scale for hopping, then we have approximately, u n « (a/k) |C n (t)| 2 . Inserting this into @, 
leads to the DNLS equation ([l]), with x = o? jk. Perhaps the most striking feature of DNLS 
is that it leads to selftrapping: for nonlinearity parameter values exceeding certain thresh- 
old (whose exact value depends on the details of the initial electron state), a portion of the 
electron remains bound to a particular site(s), while the rest propagates away in a ballistic 
manner The effect is most noticeable when completely localized initial conditions are 
used. 

The ability to partially trap an electron, suggests a possible mechanism for tuning the 
electrical conductivity. It is then clear that the robustness of the selftrapping phenomenon 
exhibited by the DNLS equation, against perturbing effects commonly found in realistic 
situations must be ascertained in order to implement possible applications. In this work we 
focus on two effects: finite inertia and quantum fluctuations. In order to keep the discussion 
simple, we will work on the smallest nontrivial DNLS unit: the nonlinear dimer Q. 

This paper is organized as follows: In section 2 we review the DNLS equation for 
the symmetric nonlinear dimer and its selftrapping phenomenology. We also consider a 
generalized nonlinear dimer. In section 3, we introduce finite oscillator inertia and examine 
its consequences on selftrapping. In section 4 we examine the fully quantum mechanical 
problem, given by the two-site Holstein model. We further specialize to the symmetric 
(section 4.1) and general (section 4.2) cases focussing on selftrapping properties. We end 
with a brief discussion of the preliminary results found so far. 
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2 DNLS dimer 



In the case of two sites, system (|T]) reduces to 

VC 2 - x \d\'d 



jlC\ r .„ , o 



^Ci-xICsl 2 ^ (4) 



For completely localized initial conditions | C\(0) | 2 = 1, an exact time-dependent solution 
for Ci(t) and C2(t) can be found|5|. The solution exhibits a selftrapping transition: the 
probability at the initial site |Ci(t)| 2 , oscillates in a periodic fashion for x < 41/, with a 
period that increases with nonlinearity. At the critical value x = ^ the period becomes 
infinite and for x > the particles selftraps. The long time-averaged probability < |Ci| 2 > 
obeys [§ 

<|Cl ' >= {l/2(l + «fc) <«> 
where ilT is the complete elliptic integral of the first kind. Another closely related system, 
the nondegenerate nonlinear dimer, considers an additional site energy difference in (H) and 
can be similarly solved Q in terms of jacobian elliptic functions. It is found that there is 
a transition line in parameter space that separates the localized from nonlocalized states. 
The difference with the previous (degenerate) case is that now the selftrapped state can be 
accessed without crossing the transition line. 

We can pose the most general nonlinear dimer by incorporating, besides a site energy 
mismatch, different nonlinearity parameters: 

dC 

= SCi+VCi-xiVifCx 

dt 

■ dC 2 Tr „ ,^ |2 



dt 



Vd-X2\C2VC 2 (6) 



The transformation d = Q exp[ (1/2) (xi — X2) \d\\ 2 ] along with the normalization condi- 
tion |Ci| 2 + \d\ 2 = 1 maps (|) into 

d~C* 

i-^ = ^C l + VC 2 - X \Ci\ 2 C l 
dt 

= Vd- X \C 2 \ 2 C 2 (7) 

where x = (V^Xxi + X2) and A = 5 — (l/2)(xi — X2)- Therefore, if we are only interested 
in site occupation probabilities, the general nonlinear dimer is equivalent to the one studied 
earlier by Tsironis^]. In Fig. 1 we show a probability density phase diagram for <P>t in 
A— x space. The transition line, which starts at x = 3.09, A = —0.383, is clearly visible. 

3 Semiclassical dimer: finite inertia effects 

As we have mentioned, the DNLS equation is obtained in the limit of a negligible oscillator 
inertia in (SI), where the oscillator adjust instantaneously to the electronic presence. We 
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now examine the robustness of the selftrapping phenomenon when finite oscillator inertia 
is taken into account. 

While in the limits of both small and large mass, analytical work is possible via perturba- 
tion theory ||], we numerically integrate the dynamical equations for the full dimer system, 
equations (||) and (|3|) with n = 1,2. The electron is placed initially completely on site 1: 
Ci(0) = 1. We explore two possible oscillator initial conditions: relaxed type, where oscilla- 
tor 1 starts from rest with a displacement determined assuming a complete oscillator relax- 
ation ui(0) = (a/k) |Ci(0)| 2 = a/k,ui(0) = 0, and undisplaced, where iti(0) = = «i(0). 
For oscillator 2, 1*2(0) = = 112(0) in both cases. The coupled equations are solved using a 
fourth-order Runge-Kutta scheme, whose precision is monitored through the conservation 
of probability (norm). 

We examine the effect of nonlinearity on the time-averaged probability at the initially 
occupied site < P > T = (1/T) J Q T \Ci{t)\ 2 dt (with T > 1/V, y/M/k), as an indicator for the 
occurrence of selftrapping. Results are shown in Fig. 2. We note that, for relaxed initial 
conditions, the finite inertia of the oscillators plays practically no role in the existence and 
location of the selftrapping transition, for at least five decades in M. There is a slight shift 
in the critical nonlinearity towards smaller values, as M is increased. On the contrary, for 
undisplaced initial conditions, the selftrapping transition is obliterated when M exceeds 0.02 
approximately. Since relaxed initial conditions are the "natural" ones when the oscillators 
have small mass, we conclude that the DNLS selftrapping transition is robust against finite 
inertia effects. Similar conclusions have been reached in the case of a number of finite 
vibrational impurities embedded in a linear chain ||. 



4 Quantum dimer 

We come now to the most delicate effect: quantum fluctuations. The DNLS equation was 
"derived" for the coupled electron-phonon system by making two approximations. The 
oscillators were first treated classically. Later, an antiadiabatic limit was taken to effectively 
eliminate the vibrational degrees of freedom, which become completely correlated (enslaved) 
with the electronic amplitude. 

In this section, we drop the previous assumptions and perform a fully quantum mechan- 
ical calculation for the rather general two-site Holstein modelQ 

H = H el + H ph + -^e-ph ( 8 ) 

with 

H el = SCld-ViClCz + Cld) (9) 

ilpb = uj (<4<ii + (4a 2 ) (10) 

#e-ph = 9iClC 1 (a\+a 1 )+g 2 ClC 2 (al + a 2 ) (11) 

where 5 is a site energy asymmetry, CJ(Cj) creates (destroys) an electron on the jth site, 

a](aj) is the usual phonon creation (destruction) operator at site j, uj is the (Einstein) 
oscillator frequency, chosen to be the same for both oscillators for simplicity and gj is the 
electron-phonon coupling at site j. The nonlinearity parameter used in DNLS corresponds 
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here to Xj = d'j/u- Previous work on a single Holstein impurity embedded in a finite chain 
have displayed significant departures from the DNLS behavior|l5|]. In that case, it was 
proved that in the antiadiabatic limit, the Hamiltonian becomes equivalent to the one of a 
static, linear impurity and <P>t never converges to the DNLS profile. 

We decouple one of the phonon fields in (|8|) using the canonical transformation 

1 ( &\ 
a = -p; ( gi ax - g 2 a 2 J 

G UJ 

b = 77 ( 92 ai + gi a 2 + ^—^) (12) 
G uj 

where G = (gf+g 2 ) 1//2 is an effective coupling. Using A = 5+{g2— g 2 ) /uj as the asymmetry 
parameter, the transformed Hamiltonian reads 

H = (A + G 2 juj)C\C x + V (C\C 2 + Cjd) + ujo)cl + G C\d (a + fl t) + ( w tfb - gj/u) (13) 

The Hamiltonian thus simplified allows for an exact numerical diagonalization scheme, 
where the eigenenergies and quantum amplitudes are computed by expanding in a phonon 



basis set[lC]. With the eigenenergies and eigenfunctions we are ready to compute dynamical 
observables. As the initial condition for the electron we choose the one that places it 
completely on one of the sites ("site 1") and focus on the probability P(t) for finding it 
there at an arbitrary time later. As before, we shall pay particular attention to the long- 
time average of P(t), <P>t= 0-/T) Jq P(t)dt, as well as <P>oo= liniT-»oo <P>t- For 
the phonon part, we use two different initial conditions: undisplaced oscillators (i.e., zero 
phonons present) and oscillators naturally "relaxed" with the electron presence (a coherent 
state) . 



4.1 Symmetric quantum dimer 

Let us consider the special situation A = 0. In this case it is possible to prove that 
<-P>oo= 1/2, independent on the initial conditions [ficfl- The original Hamiltonian can be 
separated into two diagonal blocks H± with eigenvalues {e±} and eigenvectors {|e± >}. The 
probability P(t) reads 

P(t) = (1/2) + £ C(|e± >) cos((e + - e_)t) (14) 

where the coefficients C(|e± >) depend on the initial conditions and on the eigenvectors 
|e±>. Barring any degeneracy, <P>oo= 1/2, regardless of the initial conditions. 

In the rest of this section we will use the coupling parameter g = Gj \/2 since it allows a 
better comparison with previous works. It corresponds to the coupling per site when 5 = 0. 

We are particularly interested in the adiabatic (u> <C V) and antiadiabatic {uj S> V) 
limits. For the antiadiabatic case, we obtain 

P(i)«(l/2)(l+cos(2y eff t )) (15) 

where the effective hopping term V e Q decreases with g/uj as V e g = V exp(— a(g/uj) 2 ) (with 
a = 0.434 for uj > 10). This is the well-known "polaronic narrowing" of the band, a limit 
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well understood which can be obtained by first performing a Lang-Firsov transformation 
on the Holstein Hamiltonian (||), followed by the u/V 3> 1 limit. For the adiabatic regime, 
the behavior of the electron is similar to the previous case only in the limits g/u < 1 
or g/u S> 1. In the intermediate region, P(t) displays a complex behavior due to the 
simultaneous presence of many different frequency scales pi]] 

Even though <-P>oo= 1/2 would seem to preclude any selftrapping phenomenon like 
the one observed in DNLS, it is clear from ([Tq), that the effective hopping can decrease to the 
point of "trapping" the electron for long times on the initial site. For a given observational 
time scale T, we examine < P >t- Results are shown in Figs. 3a and 3b. We note that 
the behavior of < P >x is analogous to its DNLS counterpart [12], with a "critical" value 



of g/u above which the electron tends to confine itself at the initial site, forming a "small 
polaron" , a localized electron with low mobility. 

The transition region between free electron and a small polaron behavior, can be deter- 
mined by the following criterion [fL"3j|: examination of the correlation xi,i =< ( Ao| CiCi( a 2 + 
0-2) \4>o> where \4>o> is the ground state. As g/u is increased, decreases initially from 
zero, reaches a minimum value and rises again back to zero. The maximum value of its 
slope dxi,i/d(g/u) determines the critical point {g/u) c and the maximum (minimum) value 
of its second derivative d 2 xi,i/d(g/u) 2 mark the lower (upper) boundary of the transition 
region. This is shown in Fig. 3c. Other criteria, such as variations of the ground state 
energy or the examination of xi,o =< 4>o I C\C\{a\ + a\) \4>o > give the same qualitative 
results. The features displayed in Fig. 3c seem to remain valid also for longer chains, as 
shown by Capone et al.fl3|l for N = 4 sites using exact diagonalization and by Romero et 



al. 1 14] for N = 32 sites, who relies on a variational scheme to obtain |0q > plus an energetic 
criteria. The latter group propose for the "critical" line the relation (g/u) c = 1 + (V/u) 1 ^ 2 . 
The agreement for the dimer case is good, as can be seen in Fig. 3c. According to Capone 



et al.[13], a criterion for the formation of a small polaron is the simultaneous occurrence of 




- ) > 1 and ( I > £ (16) 



that is, there must be a net phononic displacement g/u correlated with an energetic gain 
g 2 ju from polaron formation greater than the kinetic energy (~ 2V) lost due to electron 
localization. 



4.2 General quantum dimer 

The general quantum dimer is described by two parameters: the asymmetry A = 5 + (g 2 — 
g 2 )/u and the effective coupling G = ( g\ + g\ ~) x ' 2 . We focus on the values for < P > CXl 
in terms of A and G, and for the two phonon initial conditions: undisplaced and relaxed. 
Results are displayed in Fig. 4 in the form of probability density phase diagrams. 

A comparative study of < P >oo in terms of A requires to consider two scales. One of 
them is the hopping V. For small G/u values, variations in <P> 00 are determined by this 
scale. The other relevant scale is the frequency u. When the phonons are strongly coupled 
to the electron (i.e., when G/u is large enough), the effective electronic hopping increases 
when the asymmetry value is a multiple of the frequency: A ~ nu with n = 0,1, 2, ... for 
relaxed initial conditions and n = 0, ±1, ±2, ... for undisplaced initial conditions. On these 
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"resonance" channels, <P>oo= 1/2. As uj is decreased, these channels become narrower 
and begin to distort considerably the low asymmetry region (|A| < 5) giving rise to a 
probability pattern of great complexity. One of its most salient features is the presence of 
basins or low-probability "pockets" where <P>oo is significantly smaller than 1/2. These 
pockets appear for A > and can be clearly appreciated in Fig. 5, which shows a blow-up 
of the density phase diagram for the adiabatic case with relaxed initial conditions. Inside 
one of the "pockets" , the phonons conspire to push the electron towards the other site, the 
one with lower site energy. On the contrary, for the adiabatic case with undisplaced phonon 
initial conditions, no such pockets have been found and, in fact, <P>oo is always greater 
or equal to 1/2. 

Figure 5 also shows a blow-up of the probability density phase diagram for the DNLS 
case. We observe, near the DNLS critical line separating localized from non-localized states, 
a low-probability pocket, whose origin is uncertain at this stage. 

The onset of the secondary resonances (n 7^ 0) tends to occur for ever increasing values 
of G as uj decreases and, for fixed uj, as n increases. The beginning of fully consolidated 
channels can be used to define the position of a critical coupling value G c = G c (uj, A/u). 
Above this value, the average probability scales as <P>oo= ^(A/uj) where f is a smooth 
periodic function ranging from 1/2 to 1. We believe this is the small polaron region for 
the general case. A finite-time average < P >t of order unity in this region would confirm 
this conjecture. This would imply a DNLS-type behavior for < P >t- It remains to 
be determined whether the G c line approaches the DNLS selftrapping line (Fig. 1). If 
true, DNLS could be consider to appropriately describe the selftrapping dynamics of the 
quantum case, at least qualitatively, when finite-time average is considered. We are currently 
endeavoring to explain and model the resonance channels observed. They have also been 
seen for the case of a single vibrational impurity in a finite linear chain [jL3]. 

One of the authors (C.A.B) acknowledges partial support from a FONDECYT doctoral 
grant (project 2980033). 
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Captions List 



Fig. 1 : DNLS dimer: probability density phase diagram of < P >t in asymmetry- 
nonlinearity space (T = 300 V^ 1 ). The transition line begins at x = 3-09, A = —0.383. 

Fig. 2 : Semiclassical dimer: Time-averaged probability at the initially occupied site 
< P >t (vertical) versus nonlinearity parameter a 2 /k (horizontal) for different oscillator 
masses. Left column: Relaxed oscillator initial conditions. Right column: Undisplaced 
oscillator initial conditions (T = 200 V^ 1 , k = 1). 

Fig. 3 : Symmetric quantum dimer: The transition region from free electron to small 
polaron behavior is depicted by the shadowed area in (c). In (a) and (b) we show represen- 
tative <P>t plots in the adiabatic and antiadiabatic region, respectively (T = 500 V^ 1 , 
V = 1). The thick (thin) line corresponds to initially relaxed (undisplaced) oscillators. 

Fig. 4 : General quantum dimer: Probability density phase diagrams for < P >oo in 
asymmetry- "nonlinearity" space. Upper row: Adiabatic case (u> = 0.5). Lower row: An- 
tiadiabatic case (to = 10). The left (right) column refers to relaxed (undisplaced) initial 
oscillators condition. (V = 1). 

Fig. 5 : Blow-up of the probability density phase diagrams for the semiclassical nonlinear 
dimer (upper) and the general adiabatic quantum dimer (lower, uj = 0.5). 
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